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Abstract. Using a combination of the ladder operators of Piha [Rev. Mex. Fis. 
41 (1995) 913] and the parametric operators of Mielnik [J. Math. Phys. 25 (1984) 
3387] we introduce second order linear differential equations whose eigenfunctions 
are isospectral to the special functions of the mathematical physics and illustrate 
the method with several key examples. 

Resumen. Usando una combination de los operadores de escalera de Piha [Rev. 
Mex. Ffs. 41 (1995) 913] y de los operadores parametricos de Mielnik [J. Math. 
Phys. 25 (1984) 3387] introducimos operadores lineales de segundo orden con 
eigenfunciones que son formas isoespectrales de las funciones especiales de la fi'sica 
matematica y presentamos algunos ejemplos basicos. 



1 Introduction 

The use of the factorization method [T] proved to be a powerful tool for extending the class 
of exactly solvable Sturm-Liouville problems especially in quantum mechanics, where in 
the form of super symmetric quantum mechanics led to new potentials, which are isospec- 
tral to a given problem [2] . In a paper by Mielnik |3. , the usual factorization operators 

of the one dimensional harmonic oscillator Hamiltonian 

1 1 d 2 1 „ 1 

H + - = — + -x 2 + - 2 

2 2 dx 2 2 2 v ; 

have been replaced by new operators 

In order that these new operators factorize the same Hamiltonian (J2J 

bb* =aa* = H+^ (4) 

the function (3{x) should satisfy a Riccati equation of the form 

P' + P 2 = 1 + x 2 . (5) 



Using the solutions of this equation 

P(x) = x + 



7 + /e v 2 dy 
o 



(6) 



one can introduce a new Hamiltonian H, which is defined by the inverse factorization of 



b*b = H-- 
2 



with new potential functions 



2 dx 



~f + fe~y 2 dy 
o 



v(x) = — -L 

2 dx 



~f + fe~y 2 dy 
o 



and whose eigenfunctions 

i>n = b*if) n -i (n = 1,2, ...) 
are isospectral to the harmonic oscillator eigenfunctions ifi n . 



(7) 



(9) 



2 Factorization of special functions 

The great majority of differential equations appearing in Physics can be factorized by 
means of ladder operators. Therefore, we should be able to apply the procedure de- 
scribed in the previous section to the raising and lowering operators of the important 
class of Sturm-Liouville problems and get in this way isospectral second order differential 
equations. To attain this objective we proceed as follows. 

One of the possible forms of factorizing a subclass of second order differential operators 
associated to the special functions of mathematical physics was introduced by Pina [3]. 
Consider the Sturm-Liouville problem 



d 2 d 
P{x)— + Q{x)— + R n {x) 
ax z ax 



ip n (x) = 



(10) 



where P, Q and R n are functions of the variable x, and R n depend on the index n. Then, 
it is possible to construct raising and lowering operators jl] 



At 



4 + < 



a: 



P Tx + a « 



2 



that can factorize Eq.flTU]) in two ways 



A~ A + 
A + A~ 



L + K 



(12) 
(13) 



where the constant K n is the same in both factorizations. Furthermore, the functions a, 

Q d 



a n+ i, turn out to be 



A n+1 



Q d 



f 1 

^/P + c n + J —={R n+1 - R n )dx 



d f 1 

~ dx ~ Cn ~ J "7p^ n+1 ~~ Rn ^ dx 



(14) 



(15) 



<P dX J yJP 

where c n is an integration constant. From Eq.(|12Jl. one may consider the constant K n as 
the eigenvalue corresponding to the eigenfunction tjj n for the operator A~ +l A+. 
Let us now define new operators B~ +1 by 



B+ = A+ + b+ 
B- = A~ n + b- 



and demand that they can also factorize the Sturm-Liouville operator C n (x) as 

B n+ \Bn = A n+1 A+ = C n (x) + K n . 
Then, the following relationship should be fulfilled 

(b+^ + K+i^j + v^P^ + b+a~ +1 + a+b~ +1 + b+b~ +1 = . 

Therefore, the functions b£, &n+i mu st satisfy 

K+i = -K 



- - o . 

Eq.(|2*njl is a Riccati type equation, which can be easily solved to get 

b+{x) = - 



6(x) 



7 J y/P(y) V 

Here, 5(x) is defined by the indefinite integral 

{a+(y) - a~ +l (y)) 



5{x) 



P(y) 



dy 



(16) 



(17) 



(19) 
(20) 

(21) 



(22) 



and 7 is an integration constant, xq may be chosen as the point where the integrand 
vanishes. 

Similarly to Mielnik's new Hamiltonian H, we now introduce the second order differ- 
ential operator C n (x) given by 

~ db + 

C n+1 = B+B- +1 -K n = C n+l - 2y/P-^- ■ (23) 

That is, the new operator L n differs from C n by the derivative of the solution of the Riccati 
equation (|2()jl . in the same way as Mielnik's Hamiltonian H differs from the harmonic 
oscillator Hamiltonian H, as can be seen in Eq.flJJ). 
If we now define the functions 

$ n+1 = B+i> n n = 0,l,2,... (24) 

where vfj n are the eigenf unctions of C n , we can see that 

= (B+B- +1 - K n ) B^n = B+£ n i(j n = (25) 

Therefore, these ip n are the eigenfunctions for the new operator C n , and we can write the 
eigenvalue equation 

B n B n+l^n+l = + K n) = K n ljj n+1 (26) 

As can be seen from Eq. f)21jl . we have constructed a one paramenter family of operators 
£ n (x;7) which are "isospectral" to the original Sturm-Liouville operator C n (x), for the 
allowed values of the parameter 7 that produce a non-divergent function b^(x). 

As we can see, B£, B n+1 are not ladder operators as are A£, A n+X . However, one can 
easily verify that the third order operators 

Cn = B n A n _ 1 B~ (27) 

C n+1 = B n-l A n B n+l ( 28 ) 

play the role of raising and lowering operators, respectively, for the functions tp n (x). The 
use of ladder operators of order higher than two is not easy to find in the literature, but 
some work in this direction has already been reported 0. 

Since B~ +1 , and therefore C~ +1 , is not defined for n — 0, the space defined by the eigen- 
functions ifj n lack the element with n = 0, similarly to what happens in SUSY factorization 

m 

3 Examples 

Here, we proceed to find the new operators £ n (x;7) and their eigenfunctions, from the 
factorizations of the special functions of mathematical physics. 
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3.1 Hermite polynomials 

The Hermite differential equation 



d 2 H n (x) dH n (x) 



-2x- 



+ 2nH n (x) = 



dot; dec 
has raising and lowering differential operators given by 

(±- x -2 X ) Hn{x ) = -H n+1 



^-H n+1 {x) = 2{n+l)H n {x) . 
ax 



(29) 

(30) 
(31) 



In this case, the 5-integral in Eq. (j22j) is 

X 

6 = J (-2y)dy = -x 2 



and therefore 



6+ 



7 - Je~y 2 dy 
o 



(32) 



The integrand in the denominator being positive definite, one should impose the condition 
ItI > in order to have a well-defined operator. 
Now, with the use of Eq. (|2U|) one gets 



dx 



7 - Je~ y2 dy 
o 



2xe 



7 ~ Je~ y2 dy 
o 



(33) 



and therefore, the second order differential operator 

7: , \ d 2 n d Axe~ x 
C n+1 (x; 7 ) = — - 2x— + 2n + - x 



2e 



-2x z 



7 - / e v dy 
o 



7 - Je~y 2 dy 
v o 



(34) 



has parametric eigenfunctions given by 



H n+ i(x;j) = -H n+1 (x) + 



-HJx) . 



7 - Je~y 2 dy 
o 



(35) 
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3.2 Laguerre polynomials 

For the Laguerre differential equation 



x2 ^f^ + [{a + l)x ~ x2]< ^T + nxLl{x) = ' (36) 



the raising and lowering operators are 



[x^-n- L« +1 (x) = -(a + n + l)L«(x) . 



(37) 
(38) 



The 5-integral is 



and hence 



5 = J a + 2(n + 1) ~% = In x^ 2 ^ - X (39) 

^ y 



M = 5 • (40) 

7 - / y {a+2n+1) e~ydy 
o 

For x > 0, the integral in the denominator is positive definite, with maximum value 

oo 

J y (a+2n+l) e -y dy = + 2n + 2) 



and, since there is no upper limit for increasing n, we must have 7 < 0. 
The second order differential operator 



2 d 2 r , ,1 d 2b- (a + 2n + 2)lx a+2ri+1 e- a; 

J — + nx + — 5 

dx 7 -/ y^+^e-ydy 



£n+l(x;i) = x2 J^2 + (a + l)x-X 2 



9 p -2x 

r-^, < 41 > 

7 - / y a+2n+1 e-ydyj 



has as parametric eigenfunctions 



£n + i(*; 7) = (n + l)^ +1 (x) + 5 L* . (42) 

7 - / y a+2n+1 e-ydy 
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3.3 Legendre polynomials 

The Legendre differential equation 

2 i \2 d Pn{x) 



dP n (x) 



_ iyz^fL + 2x (x 2 - _ n ( n + i)(^ _ i) Pn(x ) = o 



da; 2 

has the raising and lowering operators 



d 



dx 



x 2 - 1)— + (n + l)z 



(ire 



P n (x) = (n + l)P n+ i(x) 



( ^-i) - (B+ i) 



a; 



P,H-i(a;) = -(n + l)P n (x) . 



The 5-integral is in this case 

X 

S = -2(n + 1) / —^dy = ln(l - x 2 ) 

J i-y 

(1 - x 2 ) n+1 



2\n+l 



hence 



bt 



7- /(f-?/ 2 )"^ 



(2ra + l)!! 



The integral in the denominator is positive definite, which, for a given n, has 
value 

+ 1 k/2 

J(l-x 2 ) n dx = 2 J sen 2n+1 6d6 = 2(2n) " 

-1 o 

and, therefore, I7I > 2. 

Hence, the second order differential operator 

~ d 2 d 

C n+1 (x; 7) = (1 - x 2 )— -2x(l- x 2 ) + + 1)(1 - x 2 ) + 



4(n+ 1)(1 -x 2 ) n+l 2(l-x 2 ) 2n+2 



7 - / (1 - y 2 )"^ 



7 - / (1 - 2/ 2 ) n */ 

v -1 J 



has as parametric eigenfunctions 



(1 _ x 2 ) n+1 

P n+ i{x] 7) = ~{n + l)P n +i(x) + 4 - P„(x) 

7 - / (1 - 2/ 2 ) n */ 



3.4 Chebyshev polynomials 

For the Chebyshev differential equation 

(1 " X 2 ) 2 ^M _ x{l _ X 2 } ^M + „2 (1 _ x 2 )Tn{x) = 

aar ax 
the corresponding raising and lowering operators are 

d 



[1 — x )— nx 



dx 



d 



(l-x 2 )- + (n + l)x 



dx 



T n (x) = -nT n+1 (x) . 



T n+1 (x) = (n + l)T n (x) . 



In this case, the 5-integral is 



S = _(2n + 1) f j^dy =(n+\) - 



and, hence 



bt 



(l- x 2)n+k 



7-/(l-2/ 2 ) n 2 dy 
i 



The integral in the denominator has maximum value 

1 tt/2 



j{l-x 2 ) n ^dx = 2 J sen 2n 6d6 = 



yr(2n- 1)!! 
(2n)!! 



which imposes the condition 7 > it. 

We can thus construct the second order differential operator 

C n+1 (x; 7) = (f - x 2 )^- 2 - x (f - x 2 ) L + n 2 (l - x 2 ) + 



dx 2 

2(2n+ l)x(l -x 2 ) n+ % 

i 1 

7- I(^-y 2 ) n ~ 2 dy 
-1 

with the parametric eigenfunctions 

r n+ i(x; 7) = -nT B+ i(x) + 



dx 

2(l-x 2 ) 2n+1 



-1 y 



1 - x 2 ) n+ 5 



7- / (I -y 2 )" 2 % 

-1 
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3.5 Jacobi functions 

The equation denning the Jacobi functions 



2/i fn{x) 
X (1 — X) : 

ax 



+ x(l -x)[\-(a+ l)x]^j^- + x(l - x)n(n + a)f n (x) = (57) 



has the raising and lowering operators 



/-, \ d ( n + A 

x(l — x)- [n + a) \ x — - 

v 1 dx v ; I 2n + a + 1 



Hence, the 5-integral becomes 

u — vy 



dx 



/n+lW = _ (» + l)(» + l + a-A) /t|(x) ^ 



2n + a + 1 



(59) 



= -U In x + (v — u) ln(l — x) , 



where 



u 



y(i - y) 

(n + a)(n + A) + (n + l)(n + 1 + a - A) 



2n + a + 1 



■y = 2n + a: + 1 



Therefore, we find that 



bt 



x u {l - xp~ M 



(60) 



7 _ J yu~l(l - y)v-u-l dy 

where the integral in the denominator can be written in terms of the incomplete Beta 
functions B_i(u, v — u) and B x (u, v — u). In the particular case v > u, to obtain the range 
of allowed values for the 7 parameter, we see that 



1 — 1 1 

J x u -\l - xf-'^dx = - J x u -\l - xy^dx + J x u -\l - x) v - u - l dx < 

-1 

< 2 ]x u -\l - xY~ u - l dx = 2B(u, v-u) = - T{V " U)T{U) 



T(v) 



(61) 



Therefore, we have to demand I7I > 2r ^ r ^ r( -"^ . 
The second order differential operator 



d 



d 



C n+ i(x;~/) = x (1 - x )-r-^ + x(l — x) (A - (a + l)x) — + x(l — x)n(n + a) + 



+ 



dx 2 

2{vx - u)x u {l - x) v ~ u 

7 - J y u ~ l (l - y y-v-Hy 

-1 



dx 

2x 2u {l - x) 2(v ~^ 



7- / y u -\l - y)w-Hy 
\ -1 j 



(62) 



has the parametric eigenfunctions 



7 i \ {n + a){n + A) 
^( g ^)= 2n + a + l fn+l{x) + 



x u {l-x) v ~ u 



-i 



3.6 Jacobi polynomials 

The differential equation defining the Jacobi polynomials 



dx 2 



+ (l-x 2 )[/?-a-(a + /3 + 2)x] 



+ 



+ (1 - x 2 )n(n + a + (3 + l)Pf(x) = 
has raising and lowering operators given by 



(1-x 2 )— + (n+l + a + P)[-x + 



(3 — a 



dx 



2n + 2 + a + f3 



2(n+ l)(w + 1 + a + (3) 
2n + 2 + a + (3 



(1-x 2 ) — + (n+l)(x + 
ax 



[3 — a 



2n + 2 + a + (3 



P n +\{ x ) — 



= 2(n+l + a)(n + l + (3) p 

2n + 2 + a + (3 n 1 ' ' 

In this case, we have 

/P — 92/1 1 
1 _ 2 ^ = -(g + p)ln(l + x) + -(g-p)ln(l-x) , 



where 



Hence 



P 



/3 2 -a 2 



<2 = 2n + 2 + a + /3 



2n + 2 + a + /3 

(l + x)^ q+p \l -a;)5(9-p) 
7 - J (1 + y)^^)-^! - y^-P^dy 



For the parameter 7, since q > p, we demand that 



7 > 



/ (1 + aOsfo+ri-^l - x)^^)" 1 ^ = 2"- 1 1 2 J 1 2 J . 
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From here, we construct the second order differential operator 

C n+1 (x;i) = (l-x 2 ) 2 ^ + (l-x 2 )[(3-a-(a + (3 + 2)x}l- + 

2W t , 2(qx - p)(l + x)^ q+p \l - z)3(9-p) 

+ (l-x 2 )n(n + a + /? + !) + 



7 - / (1 + y)5(<?+pM(l - y^ii-P^dy 
-l 

2(1 + x) 9+p (l - 



- - /(l + t/^+p)- 1 ^ - y )\{i-p)-^dy\ 



whose parametric eigenfunctions are 



pa/3 , x 2(n + l)(n + l + a + /3) ^ 

+ (l + s)*fr+*)(l-s)*(«-ri ^ ^ 

7 - J (1 + y)^)- 1 ^ - y^i-P^dy 
-l 

3.7 Bessel functions 

For the Bessel differential equation 

d 2 Jn(x) + 1 dJ w (x) + _ lA = o 
ofo; 2 x dx x 2 n 

the raising and lowering operators are 

(^ + l±±)j r + 1 (x) = J n (x). 



The 5-integral is found to be 



J (-1-1+1) dy = In x-^+V, 

J V v y J 



and therefore 

x -(2n+l) 2n 



n x 



7'- / y-(^+i) d y l x2n+l + * 

oo 

with 7 = 2n7 ; > 0. 



11 



From here, we can construct a second order differential operator defined by 
~ , , d 2 Id ( (n + l) 2 \ 4n + Anj(2n + l)x 2n ~ l 

whose parametric eigenfunctions are given by 

In 

Jn+l{x\ 7) = -Jn+l(x) + 2w+1 Jnfo) ■ (75) 

Hence, in the case of Bessel functions the new eigenfunctions J n+ i(x;j) are not regular 
at x — 0, except in the case 7 = 0, for which C n+1 {x]^) = and J n +i{.x',l) = J n -i( x )- 

4 Conclusion 

We elaborated here on a combination of a class of Sturm-Liouville ladder operators and 
one-parameter operators of Mielnik type that allowed us to construct isospectral Sturm- 
Liouville second-order linear differential operators with parametric eigenfunctions. Cal- 
culations are worked out for a few important cases. 
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